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ABSTRACT 

The absence of recognizable, low energy quantum gravitational effects re- 
quires that some asymptotic series expansion be wonderfully accurate, but 
the correct expansion might involve logarithms or fractional powers of New- 
ton's constant. That would explain why conventional perturbation theory 
shows uncontrollable ultraviolet divergences. We explore this possibility in 
the context of the mass of a charged, gravitating scalar. The classical limit of 
this system was solved exactly in 1960 by Arnowitt, Deser and Misner, and 
their solution does exhibit nonanalytic dependence on Newton's constant. 
We derive an exact functional integral representation for the mass of the 
quantum field theoretic system, and then develop an alternate expansion for 
it based on a correct implementation of the method of stationary phase. The 
new expansion entails adding an infinite class of new diagrams to each order 
and subtracting them from higher orders. The zeroth order term of the new 
expansion has the physical interpretation of a first quantized Klein-Gordon 
scalar which forms a bound state in the gravitational and electromagnetic 
potentials sourced by its own probability current. We show that such bound 
states exist and we obtain numerical results for their masses. 
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1 Introduction 



The problem of quantum gravity is that perturbative loop corrections to 
quantum general relativity contain ultraviolet divergences that can only be 
absorbed by adding higher derivative counterterms which would make the 
universe decay instantly [1]. The divergence problems are well known and 
ubiquitous: 

• Einstein + scalar requires a bad counterterm at one loop order [2]; 

• The same unacceptable one loop counterterm is also needed for Einstein 
+ Maxwell [3], Einstein + Dirac [1] and Einstein + Yang-Mills [5]; 

• The Einstein theory by itself requires an unacceptable counterterm at 
two loop order [6]; and 

• Although supergravity is on-shell finite at two loop order [7], and ex- 
plicit computation shows that N = 8 supergravity is on-shell finite at 
three [8] and four loop order [9], all supergravity models are expected 
to require unacceptable counterterms by seven loop order [TOj . 



Quantum gravity can of course be used as an effective field theory by 
treating the bad counterterms as perturbations and then restricting to low 
energy predictions [H] which are insensitive to them. If the Asymptotic 
Safety scenario ^12] is realized, it might even be that the escalating series of 
perturbative counterterms does not spoil predictivity at energies below the 
Planck scale. However, neither approach provides a fundamental resolution. 

The problem arises from the tension between four facts [T]: 

• Continuum Field Theories possess an infinite number of modes; 

• Quantum Mechanics requires each mode to have a minimum amount 
of energy; 

• General Relativity stipulates that stress-energy is the source of gravi- 
tation; and 

• Perturbation Theory simply adds up the contribution from each mode 
at lowest order. 
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One or more of these principles must be sacrificed, and a little thought sug- 
gests focussing on the last two. There does not seem any way of disputing 
the experimental confirmation of quantum mechanics in the matter sector 
which is responsible for the lowest order divergences of quantum general rel- 
ativity. And inflationary cosmology makes nonsense of any attempt to invoke 
a nonzero physical cutoff length. Inflation predicts that the universe has ex- 
panded by the staggering factor of at least 10^^ [13], so if the physical cutoff 
is at the Planck length today then it must have been about 10~^° m dur- 
ing primordial inflation. But fossilized quantum gravitational effects from 
primordial inflation have been measured with a fractional strength of about 
10^^° [11], which is inconsistent with so small a physical cutoff length. 

Superstring theory can be viewed as an attempt to preserve the validity of 
perturbation theory by sacrificing general relativity. We wish here to investi- 
gate the alternative: that the problems of quantum general relativity derive 
from using conventional perturbation theory. We disavow any intention of 
seeking the exact solution. There is so far no example of an interacting quan- 
tum field theory in D = 3 + 1 dimensions which can be solved exactly, and 
all experience with classical field theory suggests that general relativity is an 
unlikely candidate to be the first one. What interests us instead is the possi- 
bility that quantum general relativity has a perfectly finite asymptotic series 
expansion which is simply not given by conventional perturbation theory. 

The conventional result for the expectation value of a quantum gravity 
observable O with characteristic length R is assumed to take the form. 



O 



r]) = (TreeOrder)|l + f:a,(J^)1, (1) 



where G is Newton's constant, % is Planck's constant and c is the speed of 
light. Support for this form is adduced from the fact that quantum grav- 
ity has no observable effects at low energies. Even for the smallest dis- 
tances ever probed, R ~ 10~^^ m, the loop counting parameter is minuscule, 
TiG/{c^R^) ~ 10~^^. But the same thing would be true of a series that 
incorporates fractional powers or logarithms such as. 



O 



hC 



i=l m=0 ^ 



hG 



O) = (T.ee Orie.)] 1 + E ZMS) H^T \ ■ 



If the actual asymptotic expansion of quantum gravity were to take the form 
(I2l) then loop effects at i? ~ 10~^^ m would still be suppressed by unobserv- 
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ably small powers of the parameter, 



/ TiG \ / TiG 



~ -10-^" . (3) 

R~10-19 m 



However, trying to force the putative series (|2]) into the assumed form ([T]) 
would result in logarithmically divergent coefficients a^, which is exactly what 
explicit computations reveal. 

The incorporation of such nonanalytic terms into an asymptotic expan- 
sion occurs even for very simple physical systems. Consider the canonical 
partition function for a non-interacting, highly relativistic particle of mass 
m in a three dimensional volume V at temperature UbT = 1/13, 



y /"OO r / -1 

— — ^ / dpp^exp —pJp^c^+rri^d^ + Pmc^ , (4) 

2TT^h JO L V J 
T/ roo 

^ o , / dK (K+mc^) yjK'^ + 2Kmc^ exp{-f3K) . (5) 

2TT^h Jo 



When the rest mass energy is small compared to the thermal energy it ought 
to make sense to expand in the small parameter x = (imc^. But straightfor- 
ward perturbation theory fails. 



From expression ([7]) it seems as though the term of order x^ vanishes, and 
that the higher terms have increasingly divergent coefficients with oscillating 
signs. In fact the x^ term is non-zero, and the apparent divergences merely 
signal contamination with logarithms, 

Just as we suspect is the case for quantum gravity, the partition function 
has an excellent expansion for small x = Pmc^; the terms after order are 
indeed smaller than x^, they just are not as small as one naively thinks. 

Rather than attempting to develop a new expansion for an arbitrary 
quantum gravity observable, we restrict attention here to the self-energy of a 
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charged, gravitating particle. An exact result for the ^ = limit of this sys- 
tem was obtained in 1960 by Arnowitt, Deser and Misner (ADM) [T5]. Their 
work provides strong support both for the possibility that negative gravita- 
tional interaction energy cancels divergences, and for the possibility that the 
correct asymptotic expansion involves nonanalytic dependence on Newton's 
constant. We review this evidence in section 2. In section 3 we propose an 
alternate expansion scheme for the self-energy of a quantum field-theoretic 
particle. How the new expansion reshuffles the diagrams of conventional per- 
turbation theory is worked out in section 4. We discuss the 0th order term 
of the new expansion in section 5. Our conclusions comprise section 6. 



2 The ADM Computation 

Arnowitt, Deser and Misner showed that perturbation theory breaks down in 
computing the self-energy of a classical, charged, gravitating point particle 
[15j. It is simplest to model the particle as a stationary spherical shell of 
radius i?, charge e and bare mass Mq. In Newtonian gravity its rest mass 
energy would be, 

= + - ^ • 

It turns out that all the effects of general relativity are accounted for by 
assuming it is the full mass Mpi which gravitates, rather than Mq [TS] , 



e2 GMl m 2G / _ e2 ^ Rc^ , ^ 

The perturbative result is obtained by expanding the square root, 

P " 87reoi?c2 ^' \ Rcy \ SireoRc^ J ' ^ ^ 

and shows the oscillating series of increasingly singular terms characteristic 
of the previous examples. The alternating sign derives from the fact that 
gravity is attractive. The positive divergence of order e'^/R evokes a negative 
divergence or order Ge^/R^, which results in a positive divergence of order 
G'^e^/R^, and so on. The reason these terms are increasingly singular is that 
the gravitational response to an effect at one order is delayed to a higher 
order in perturbation theory. 
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The correct result is obtained by taking R to zero before expanding in 
the couphng constants and G, 



Like the example of Section [T] it is finite but not analytic in the coupling 
constants and G. Unlike this example, it diverges for small G. This is 
because gravity has regulated the linear self-energy divergence which results 
for a non-gravitating charged particle. 

One can understand the process from the fact that gravity has a built-in 
tendency to oppose divergences. A charge shell does not want to contract in 
pure electromagnetism; the act of compressing it calls forth a huge energy 
density concentrated in the nearby electric field. Gravity, on the other hand, 
tends to make things collapse, especially large concentrations of energy den- 
sity. The dynamical signature of this tendency is the large negative energy 
density concentrated in the Newtonian gravitational potential. In the limit 
the two effects balance and a finite total mass results. 

Said this way, there seems no reason why gravitational interactions should 
not act to cancel divergences in quantum field theory [TB]. It is especially 
significant, in this context, that the divergences of some quantum field the- 
ories — such as QED — are weaker than the linear ones which ADM have 
shown that classical gravity regulates. The frustrating thing is that one can- 
not hope to see the cancellation perturbatively. In perturbation theory the 
gravitational response to an effect at any order must be delayed to a higher 
order. This is why the perturbative result flTT]) consists of an oscillating 
series of ever higher divergences. What is needed is an approximation tech- 
nique in which gravity knows what is happening in the gauge sector so the 
gravitational response can keep pace at the same order. 

A final point of interest is that any finite bare mass drops out of the exact 
result f ll2p in the limit i? — )■ 0. This makes for an interesting contrast with 
the usual program of renormalization. Without gravity one would pick the 
desired physical mass, Mphys, and then adjust the bare mass to be whatever 
divergent quantity was necessary to give it. 




(12) 



Mo = Mphy, - 



(13) 
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Of course the same procedure would work with gravity as well, 



Mr 




(14) 



— 



phys 



87reoi?c2 2eoRc^ 



The difference with gravity is that we have an alternative: keep Mq finite 
and let the dynamical cancellation of divergences produce a unique result for 
the physical mass. The ADM mechanism is in fact the classical realization of 
the old dream of computing a particle's mass from its self- interactions |17j . 

3 A New Expansion for Particle Masses 

The purpose of this section is to explain the new expansion we propose for 
particle masses. For simplicity we work in the context of a charged and 
gravitating scalar field, although the same technique applies to fermions and 
to Yang-Mills force fields. The Lagrangian is the sum of those for general 
relativity, electrodynamics and a charged scalar. 



£sc[^?,A0,0l = -{D^4>yD,<Pg^^^g-Ml(l>*<P^g. (17) 



Here and henceforth g^jj{t,x) stands for the metric field, with inverse g^'^ 
and determinant g; A^{t,x) denotes the electromagnetic vector potential 
with field strength F^j, = d^j_Ay — d^A^; and (p{t,x) is the complex scalar 
field. The covariant derivative operator is = + ieA^. The alert reader 
will note that the scalar Lagrangian lacks the quartic self-interaction that 
would be required for perturbative renormalizability in flat space. Because 
the charged scalar is anyway not perturbatively renormalizable once gravity 
has been included, there does not seem to be any point to including this term 
for a first investigation of nonperturbative renormalizability. 

We employ the usual units of particle physics in which h = 1 = c, so 
that time and space have the dimensions of inverse mass, the charge e is a 
pure number, the Newton constant G is an inverse mass-squared, and the 
bare mass Mq is a mass. We shall also sometimes distinguish time and space 
arguments — as in (f){t, x) — and sometimes lump them together into a single 
spacetime coordinate x'^ = {t,x) — as in (f){x). 




(15) 



IQttG 




(16) 
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Our attitude is that the physical mass M of single scalar states is some 
complicated function of the bare parameters, e, G and Mq. Our first goal is 
to derive a formal expression that would give M, assuming we had infinite 
computational ability. We then develop an alternative to the conventional 
perturbative expansion for evaluating this formal expression. 

3.1 Functional Integral Expression for the Mass 

If all interactions were turned off it would be simple to express the free scalar 
field in terms of the operators and a{k) which create and annihilate one 
particle states with wave number /c, 

0free(t,x) = |^ / -i= »*+^^^-a(^) + -i=e-«*-^^-V(fc)| . (18) 



Here ujq = \J k"^ + Mq is the free energy. We can invert relation (fT8|) to solve 
for the annihilation operator using the Wronskian = — 

a{k) = ^e'^°'WohUt,k) . (19) 



Here and henceforth, a tilde over some function denotes its spatial Fourier 
transform, 

f{t,k) = Jd^xe-'''-^f{t,x) . (20) 

In the presence of interactions it is no longer possible to give explicit 
relations such as (fT9|l for the operators which create and destroy exact 1- 
particle states. However, if we temporarily regulate infrared divergences and 
agree to understand operator relations in the weak sense then it is possible to 
write the operators which annihilate outgoing particles and create incoming 
ones as simple limits |18j . 



a-*(fc) = hm -== f^,, 0(t+, A;) , (21) 



(a'°(fc))^ = lim ^^t,_r(t_,fc). (22) 

Here Z is the field strength renormalization (defined as the amplitude for the 
field to create a 1-particle state) and u is the full energy. 



VWTIP . (23) 
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At this stage we do not know the physical mass M; it is some function of the 
bare parameters of the theory. 

Now consider single particle states whose wave functions in the infinite 
past and future are V't(^)' respectively. We can employ relations f l2Tll22|) to 
derive an expression for the inner product between these states, 



(27r)3 2Zuj 



lim e''^'+Wt 



lim e 



-iujt- 



0(t+,f)0*(t_,O)fi'"). (24) 



One way of computing the physical mass M would be to adjust it to the 
precise value for which expression reduces to, 



(27r)3 2Zuj 



(25) 



This agrees with the usual definition of the mass as the pole of the propagator, 
but it is problematic owing to infrared divergences. 

A more direct way of computing the mass is to focus on the second line 
of ( IMl) which we can express as the exponent of —i times some complex 
function C,{t+, t^,k), 



0(t+,f)r(t_,o) 



(26) 



This function C,{t+,t^, k) includes many things — the field strength renor- 
malization, finite time correlation effects from multiparticle states, and so on 
— but only the single particle energy grows linearly with the time interval. 
Dividing by the time interval and then taking it to infinity gives this energy. 



lim 

t±— )-±oo 



^{t+,t.,k) 



Vk^ + M2 . 



(27) 



Setting k = gives the physical mass we are seeking. This way of computing 
M avoids the problems of infrared divergences which affect the field strength 
renormalization but not the mass. 

It is straightforward to write (l26i) as a functional integral. 



d X e 



-ik-x 



[dg][dA][d(l)][d(l)*] 0(t+,x)0*(t_,O) e^^GR[9]+iSEM[</,A]+i5sc[9,A,<A,r] _ (28) 
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We have subsumed the details of gravitational and electromagnetic gauge 
fixing into the measure factors [dg\ and [dA\. Note that the various action 
integrals in expression fl28|l go from time t_ to t+, 



SgM = dtjSxCoMit.x) , (29) 
S^mIq^A] = j'jdtjd\C^M[9.A]{t,x), (30) 
^sc[^?,A0,0l = jytJd^xCsc[9,A,(l),(l)*]{t,x) . (31) 

3.2 Eliminating the Matter Field 

Expression f l28|) is only formal because there is not yet any way of defining 
it or evaluating it. One would usually resort to perturbation theory at this 
point. We shall instead integrate out the scalar field, 

g -i5{t+,t_,fe) _ /"(^3^ g-ifc x 

[dg][dA] 2A[^,A](t+,f;t__,0)e*^««[^^]+^^^^[^''^]+^^sc[9,A] _ (33) 

The two new quantities this produces are the scalar propagator iA[g, A] {x; x') 
and the scalar effective action Fscifi', Both can be defined using the scalar 
kinetic operator in the presence of an arbitrary metric and vector potential, 

V[g, A] = D^^g^'^D, - . (33) 

In rough terms, the scalar propagator is i times the functional inverse of 
V[g,A], subject to Feynman boundary conditions, while the effective action 
is i times the logarithm of its determinant, 

V[g, A] X iA[g, A]{x; x') = i6\x-x') , (34) 
Tsc[9,A] =t\n[det[V[g,A]-ze]) . (35) 

It will facilitate subsequent work to be more precise about the definition of 
of the scalar propagator iA[g, A]{x; x'). Consider a general "mode function" 
u[g, A](t,x] A) which obeys the homogeneous equation, 

V[g,A]xu[g,A]{x;X) = 0. (36) 
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Here "A" is a (possibly continuous) index which labels the solution; it would 
be the wave vector k for flat space and zero potential. We additionally require 
that the set of all solutions obey the canonical normalization condition, 

-I [ d'xJ-g{x) u[g, A]{x- \)(% - t*)u*[g, A]{x- k) = . (37) 



't=const 

In terms of these mode functions the propagator is, 



iA[g,A]{x;x') = Y: 



eit-t')u[g,A]ix;X)u*[g,A]ix';X) 

+9{t'-t)u*[g,A]{x;XHg,A]{x'; A) 



(38) 



Note that only the first theta function contributes in the limit we require, 
Jim iA[g,A](t+,x;t_,6) = J2u[g, A]{t^, x; \)u*[g, A\{t_,6; \) . (39) 



Our expression for the physical mass can therefore be written as. 



M = lim 



In 



[dg][dA] 



xu[g, A]{t+, X- A) u*[g, A](t_, 0; A) e 



3.3 Stationary Phase Expansion 



(40) 



The expression poj) we have derived for the physical mass is exact, but 
formal because no one knows how to evaluate the functional integral. That 
would ordinarily be done by resorting to conventional perturbation theory. 
This could be viewed as an application of a simplified variant of the method 
of stationary phase to the original functional integral fl25]) . Because our 
modified expansion involves undoing some of the simplifications we digress 
to review the technique. 

Recall that the method of stationary phase gives an asymptotic expansion 
for integrals of the form, 

/ = Idze'^^^^ . (41) 
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The technique is to first find the stationary point zq (which we assume to be 
unique) such that f'{zo) = 0. One then expands f{z) around zq, 



1 oo f{n)f \ 

f{zo) + -f"{z,){z-z,f + {z-z,r 
h + \n{z-z,f + l^f{z-Zo). 



The next step is shifting to the variable C, = z — zq and expanding e*^-^, 

oo 



J zT^n ml L 



m=0 



(42) 
(43) 

(44) 



At this stage the result is still exact, but generally no simpler to evaluate 
than the original form. What gives a computable series is the final step 
of interchanging integration and summation. It is at this point that the 
expansion ceases to be exact and becomes only asymptotic, 



oo jm 



E 

m=0 



ml 



64' 



A/(C) 



(45) 
(46) 



Conventional perturbation theory is a simplified form of this technique, 
applied to the functional integral (128|) . The two simplifications are: 

• The multiphcative factor of 0(t+, af)0*(t_, 0) is not included in the ex- 
ponential, along with the action; and 

• The stationary field configuration is assumed to be fiat space with no 
charge fields. 



A^{x) = 



(x) = . 



(47) 



These two simplification make conventional perturbation theory much sim- 
pler and more generally applicable than a strict application of the method 
of stationary phase because they remove any dependence of the propagators 
and vertices on the operator whose expectation value is being computed — in 
this case 0(t+, x)0*(t_, 0). However, simplicity is not always desirable, nor is 
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it always physically correct. In this case, those simplifications preordain that 
the eventual asymptotic expansion for M can contain only integer powers 
of G and e^. The same two simplifications also imply that the gravitational 
response to an £ loop effect in the electromagnetic sector must be delayed 
until £ + 1 loop order. 

Our modified expansion is defined by making three changes. The first 
is to integrate out the matter fields and start from the functional integral 
f HOj) . That is not so important. The important change is the second one: 
we include the factor u[g, A\{t^,x]X) x 0; A) in the exponential, 

along with the gravitational and electromagnetic actions [19]. This makes 
the eventual series vastly more comphcated but it does three desirable things: 

• It allows the apparatus of perturbation theory — propagators and ver- 
tices — to depend on the operator whose expectation value is being 
computed; 

• It permits the eventual expansion to involve complicated, nonanalytic 
functions of and G; and 

• It allows the gravitational response to an electromagnetic effect at some 
order to occur at the same perturbative order. 

The final change is that we drop the scalar effective action Fscfs', ^] from 
how we define the expansion. That is, it plays no role in determining the 
stationary point, the "classical" action, the propagators or the vertices; we 
regard e*'"^*^ as a multiplicative factor, like (fxp* was in the conventional ex- 
pansion. Although this can be done consistently because Fgcis',^] is gauge 
invariant, there is no physical justification for it. The scalar effective ac- 
tion incorporates one loop vacuum polarization and its gravitational analog, 
which may be important effects. Our justification is just the pragmatic one 
that including rscif?, ^] leads to a more complicated set of equations for the 
stationary point than we presently understand how to solve. 

3.4 Physical Interpretation 

The 0th order term in the new expansion can be interpreted as the energy of 
a first-quantized Klein-Gordon particle which moves in the electromagnetic 
and gravitational fields that are sourced by its probability current pQ]. To see 
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this, note that the quantum mechanical propagator for such a Klein-Gordon 
particle in fixed background fields A^{x) and g^u{x) is, 

V[g,A]{x-x') = Y.u[g,A]{x-\)u*[g,A]{x'-X) . (48) 

A 

That means we can evolve the first-quantized wave function from t' to t by 
taking the inner product with V[g, Al\{x; x') according to relation fl57|) . 

/ Sx'^J-g{x')g''\x')V[g,A]{x-x')W,u[g,A]{x'-\)=u[g,A]{x-\) ■ (49) 

J t'=const 

In expression (HOi) we are going from a delta function il)-{x') = (5^(x') to 
a zero momentum plane wave ip+{x) = 1. Of course the stationary field 
configurations for the metric and the vector potential are just those sourced 
by the quantum mechanical particle itself, hence the stated interpretation of 
the 0th order term. 

Our ability to evaluate the 0th order term depends upon whether or not 
the first-quantized particle can form a bound state in its own potentials. If 
not then we are left with a complicated scattering problem which seems to 
be intractable. However, many simplifications are possible if a bound state 
forms. First, one can forget about the continuum solutions; the result for M 
in expression ( HOl) will derive entirely from the bound state with the largest 
overlap between the two asymptotic wavefunctions 'ip±. Second, one can 
specialize the wave function to take the form, 

u[g,A]{t,x) = e'^'''F{r) (50) 

Third, in solving for the stationary potentials one need only consider a class 
of metrics and vector potentials which is broad enough to include the eventual 
solution. For scalar QED this reduces the potentials from nine functions of 
spacetime (after gauge fixing) to only three functions of a single variable, 

g^^{x)dx^'dx'' = -B{r)dt^ + A{r)dr^ + r^dn'^ , (51) 
A^{x)dx'' = -$(r)rft. (52) 

This simplification might also be relevant to the problem of going beyond 0th 
order. Although it is impossible to work out propagators for general back- 
ground fields, it is sometimes possible to derive the propagator for potentials 
which depend upon only a single variable [21]. 
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A final simplification is the possibility of deriving a variational formula 
for the bound state energy. Even if the functions F{r), A{r), B{r) and $(r) 
could not be determined exactly, this technique would allow us to explore 
reasonable approximations for them. We might even optimize the free pa- 
rameters in the largest class of solutions for which the field propagators can 
be worked out. A variational technique would also give an upper bound on 
the bound state energy. 

The reason deriving a variational technique is only a possibility is that 
this system involves gravitation. It is well known that not all of the gravi- 
tational potentials contribute positive energy |22]. Of course it is precisely 
the negative energy potentials that provide the possibility for canceling self- 
energy divergences! These negative energy potentials do not lead to any 
instability because they are completely constrained variables; that is, they 
are determined in terms of the other variables. However, it is the constraint 
equations which enforces these relations. The Hamiltonian is not bounded 
below before imposing these constraint equations. So being able to derive a 
variational formalism depends upon identifying the negative energy poten- 
tials and solving the constraint equations for them. 

4 Rearrangement of Perturbation Theory 

The purpose of this section is to show where the old diagrams end up in the 
new expansion. The point of doing this is not to perform the computation; 
there are more efficient techniques which will be developed in the next section. 
The point is rather to see that the new expansion offers gravity the chance 
of "keeping up" with what is happening in other sectors. We will show 
that all the usual diagrams are present, and at the same "loop" order as 
in conventional perturbation theory. However, they are joined with a vast 
class of new diagrams which are added to one order and subtracted from 
another. To see all this we resort to a very simple, zero dimensional model 
of the functional integral in expression fl28|) . We first work out what the 
conventional expansion gives, then derive the zeroth and first order results for 
the new expansion. The section closes with a discussion of the implications 
for the new expansion. 
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4.1 Zero Dimensional Model 



We can model the functional integral f l28|) as an ordinary integration over two 
A^- vectors: Xj, representing the metric and gauge fields; and y"' representing 
the scalar. A simple model for the action is, 



S{x, y) ^ \k,x,x^ + \B^'y^y' + -^Cfx^y^^y" . 



(53) 



Of course scalar QED includes interactions between two photons and two 
scalars, and the action of general relativity contains self-interactions of the 
metric, as well as infinite order interactions of the metric with the scalar and 
the vector potential. However, our model action (13^ will suffice for the pur- 
pose of understanding how conventional perturbation theory is reorganized. 
The normalization integral and the 2-point function of our model are. 



J 

jah 



J-^ Jd^x Jd^yy^y^e'^^^^'y^ . 



(54) 
(55) 



It will simplify the notation if we use the script letters Aij and B"''' to denote 
the matrix inverses of Aij and B"''', 



A- A -1 



Hk 



ab rybc 



B^'B 



5" 



(56) 



When multiplied by i these inverses are the "propagators" of this one dimen- 
sional quantum field theory. The expansions of the normalization integral 
and the 2-point function are. 



J 



(27r)^e 



1+^iCfiCfiAij iB'"'iB'"'+2iB^HB"' +0(C^) \. (57) 



2a6 • toed 



ac- t2bd 



det(A)det(5) I 8 



jab ^ iB'''' + UB''HCfiAijiCfiB'^iB'^'' + iB''HCfiAjiB'^HCfiB^^ + OiC^^^ . 

(58) 

Fig. [1] gives the diagrams associated with the expansion fl58l) . 

Following subsection 3.2, we integrate out the "matter fields" y"', 



B+CkXk 



ab 



(27r)-e 



JV JVir • 
4 * 



det{B + CiXi) 



(59) 
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Figure 1: Usual expansion for the scalar propagator. Solid lines are scalars, 
wavy lines represent photons and gravitons. 



At this stage it is useful to extract the lowest order contribution from the 
determinant, 



det{B + CiXe) = det{B) x det{I-iB ■ iCiXi) 



(60) 



It is also useful to contract into "asymptotic state wavefunctions" — the 
analogues of 4'±- These are the A^- vectors m", analogous to ■?/'+, and v^, 
analogous to We shall also include a phase so that u°'iB°'^v^ = 1. Then 
the contraction of these two vectors into the field dependent propagator, the 
analogue of iA[g, A]{x; x'), can be written as, 



u X 



1 



B + CkXk 



ab 



X V 



U" X 



I—iB ■ iCkXk 



X iB"^ X V' 



(61) 



= 1 + n'^ X ^ 

n=l 

= 1 + iA(x) . 



{iB ■ iCkXk) 



xiB^'xv', (62) 



(63) 



Except for the phase, the quantity whose logarithm we wish to take is. 



det(A) 



JV JVir • 
4 * 



(27r)-e- 

Here the exponent £{x) and the "matter effective action" T{x) are. 



(64) 



Tix) = - In 
^ ' 2 



dei{I-iB-iCkXk) ■ (65) 
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+ 



+ 



u - 
u - 



+ 



Figure 2: Expansion for the solution Xi to equation ( 1671) . Solid lines are 
y"" propagators (scalars), wavy lines represent Xi propagators (photons and 
gravitons) . 

4.2 0th Order in the Model 

We want to do the stationary phase expansion on expression (IMI) . but treat- 
ing r(x) as higher order. Hence the zeroth order term in the expansion is, 



l+iA{X) 



g 2 ^ij-^i-^j 



Jo = e 

where the "stationary field configuration" Xi is found by solving, 

■ aiA(x) 

^' . 

x=X 



d£{x) 



(66) 



(67) 



■' ' l+tA{X) 

Despite the simplicity of our model, equation ( 167|) is too difficult to solve 
exactly for general Aij, B"^^ and Cf^. However, we are only interested in a 
perturbative solution — in powers of the "interaction" C^^ — and that is 
simple enough to generate by iteration. 



• . diA(X) 
^ ''•^V dxj 



l+iA(X) 



(68) 



'iBiCjiB 


V — iAijU 


iBiCjiB 


V X u 


iBiCkiB 


viAkiU 


iBiCiiB 



u 



iB(iCjiBiCk+iCkiBiCj)iB\viAMu[iBiCe,iB\v + 0{C^) . (69) 

Fig. [2] gives a diagrammatic representation of the expansion. 

In evaluating Jq it is useful to first formally express things in terms of 
zA(X) and 



iA(X) 

diA{x) 



dxi 



u 



u 



iB-iCkXk-iB 
iB-iCkXk-iB 



V + u 
v + 2u 
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iB-iCkXk-iB-iCtXriB\v + ■ ■ ■ , (70) 
iB-iCkXk-iB-iCiXriB]v + ■ ■ ■ . (71) 



V I u 
_3 \ 



Figure 3: Expansion for the 0th order term Jq. Sohd hnes are y"" propagators 
(scalars), wavy hnes represent Xj propagators (photons and gravitons). 



Using expressions ( 166|) and (!67|) we can write j7o as, 

Jo = 1 + 2A(X)-- ^ ^ . , j,. +■ 



1 + 



2 



iB-iCkXk-iB 



8[l+zA(X)] 



f H — M 



iB-iCkXk-iB 



V + 



Now just substitute fl69|) in (1731) to obtain 



Jo = 1 + 2« 
-- u 



iB-ia-iB 



iB-iCriB 



V iAij u 



V iAij u 



iB-iCj-iB 



iB-iCj-iB 



V \ + u 



iB-iCi-iB 



X iAij u 



iB-iCj-iB-iCk-iB 



viAkfU 



iB-iCriB 



Fig. I3] gives a diagrammatic expansion of Jo- 



(72) 
(73) 



V + 0{C^) . (74) 



4.3 1st Order in the Model 

One loop effects come from the determinantal correction to the Method of 
Stationary Phase, 



A. 



d^£{X) 



Jo+i = e^^(^) X (^det 
The matrix whose determinant we must compute is 



(75) 



A 



dxjdxk 



Sik + 



diAjX) d'iA(X) -. 



l+iA{X) 
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= Sik + M, 



ik 



(76) 



We can evaluate the determinant in f l75p by first expanding in powers of tlie 



matrix 



det(/+M) 



g-iTr[ln(/+M)] 



(77) 



= 1 - ^Tr[M] + i (Tr [M]) ' + ^Tr [M^] + 0(M3) . (78) 
2 o ^ '4 

Tlie determinantal contribution to (!75|) is sufficiently complex that it is 
worthwhile to present results for each term of (1781) . As before, it is best to 
express these contributions in terms of the full solution Xj, 



-Tr[M] 



iB-iChXh-iB 



V + iAijU 



iB-zQ-iB-iCj-iB 



I 

+-{ u 
2\ 



iB-iCuXu-iB 



V \ — u 



iB-iCuXu-iB 



V lAij u 



iB-iC-iB-iCyiB 



—u 



iB-iCkXk-iB-iCeXriB 



V + iAijU 



iB\ iCi-iB-iCj-iB-iCkXk 



+iCi ■ iB ■ iCkXk ■ iB ■ iCj+iCkXk ■ iB ■ iQ ■ iB ■ iC^ ) iB 



V + 0{C^) , (79) 



;(Tr[M]) 



2 1 



U 



iB-iCkXk-iB 



V --U 



iB-iCkXk-iB 



V iAi 



xu 



iB-iCi-iB-iCyiB 



V + - \ iAij u 



iB-iCi-iB-iCj-iB 



v] +0(C^), (80) 



-Tr[M^] = -U 



iB-iChXk-iB 



V \ — u 



iB-iCkXk-iB-iC^Xi-iB v+-iAijiAki 

J 2 



iB-iCj-iB-iCk-iB 



vu 



tB[iCe-iB-tCi+tQ-iB-tCe)tB\v + 0{C^) . (81) 
The first order contribution to J in our new expansion is, 

d^8{X) 



Ji = e^^(^)|(^det 
1 



A 



-u 



iB-iCriB 



dxjdxk 

iB-lCyiB 



V iAij u 



1 



V + iAij u 



iB-iQiBiCj-iB 



(82) 



iB-iQ-iB 



V iAij u 



iB-iCj-iB 



f — 3m 



iB-iCi-iB 



viAi 
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- V u 



- 'U u 



Figure 4: Expansion for the 1st order term jTi- Solid lines are propagators 
(scalars), wavy lines represent Xj propagators (photons and gravitons). 



iB-iCj-iB-iCk-iB 



V iAkf u 



iB-iCriB 



V — u 



iB-iCriB 



xiAijU iB-iCj-iB viAkeu 



iB-iCi.-iB-iCriB 



V + u 



iB-CriB 



V iAij iAki 



xu 



iB[iCj-iB-iCk-iB-iCt+iCk-iB-iCj-iB-iCi+iCk-iB-iCriB-iCAiB 



+-iAijU 



iB-iCj-iB-iCk-iB 



viAkfU 



iB{iCriB-iCi+iCi-iB- id) iB 



+- iAii u 



iB-iCi-iB-iCyiB 



V + 0(C6) . (83) 



Fig. m gives a diagrammatic representation of the expansion. 



4.4 Implications for Our Expansion 

Let us compare the new expansion J = i7o + i7i + " ■ " with the old one, 
J' = Jo + Jl + ■ ■ ■• The first two terms of the new expansion are expressions 
(IMj) and flS^ . In contrast, the first two terms of the old expansion are. 



iB-iCi-iB-iCjiB 



(84) 



Jo = 1 , Jl = u 
It is useful to define the difference between new and old at order £, 

^J, = J,-J,. (85) 
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We can read off ATo from expression (1^ . 
1 



A7o 



-u 



-- u 



iB-iCi-iB 
iB-iCriB 



V lAij u 



iB-iCylB 



V iAij u 



iB-iCyiB 



V u 



xu 



iB-iCj-iB-iCk-iB 



viAkiu 



iB-iCriB 



iB-iCfiB 



viAr 



V + O(C^) . (86) 



Note that eacli contribution to ATo is a "tree diagram"; the expansion in 
powers of the interaction C could also be viewed as an expansion in numbers 
of external lines. From expression ([HSD we see that /\Ji possesses both "tree" 
and "one loop" contributions, 

1 



-u 



+ - 



7 



u 



iB-iC-iB 
iB-iCriB 



V iAij u 



V iAij u 



iB-iCyiB 



iB-iCj-iB 



3u 



iB-iCi-iB 



V iA: 



+u 



-u 



xu 

iB-iCriB 



iB-iCj-iB-iCk-iB 



viAij u 



iB-iCj-iB 



viAkiU 
viAkeu 



iB-iCe-iB 



V + 0{C^) , (87) 



iB-iCk-iB-iCg-iB 



iB-CriB 



V iAij iAki u 



iBHCj-iB-iCk-iB-iCi+iCk-iB-iCj-iB-iCi 



+iCk-iB-iCriB-iCMB 



1 

V H — lAij u 



iB-iCi-iB-iCu-iB 



viAkiu iB{iCi 



xiB-iCi + iCi-iB-iCe]iB v + -liAij u 



iB-iCi-iB-iCj-iB 



V + 0(C6) . (88) 



Note that the order contributions to ZX7o and Aj7i° cancel. They do 
not cancel at order C"^, but the residual terms are canceled by /\J2- So 
the new expansion does not move the old diagrams from one order to another, 
rather it adds a new class of diagrams — with more external lines — to one 
order and subtracts them from higher orders. What we have learned can be 
summarized by the following observations: 

• Each term in the new expansion can be written as an infinite series 
which begins at order C^^; 

• The lowest order term in this series expansion of each JT^ is the old 
£-loop term J^; 
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• The new terms in Ji involve diagrams with up to and including t loops; 

• The new terms at each order have more powers of the external state 
factors u"" and f and 

• The sum of all the new terms is zero. 

Although there is no guarantee that the new expansion is free of ultra- 
violet divergences, it does possess a number of desirable features. Note first 
that ultraviolet divergences are no worse in the new expansion than the old 
one because none of the new diagrams in Ji have more than i loops. Our 
simple model makes no distinction between photons and gravitons but the 
quantum field theoretic expansion will of course involve both particles. The 
new diagrams — with more powers of the coupling constants and G, but 
no more loops — are one way gravity at order I in the new expansion can 
know about an order I divergence from the gauge sector. So the new ex- 
pansion breaks the conundrum of conventional perturbation theory that the 
gravitational response to a problem at one order cannot come until the next 
order. Note finally, that the infinite sums of new diagrams at any fixed order 
of the new expansion allow the possibility of getting nonanalytic dependence 
upon G and e^. 

5 The 0th Order Term 

The purpose of this section is the evaluate the 0th order term in the new 
expansion of expression (HUj) for the scalar mass. Recall that this means 
solving the problem of a first quantized Klein- Gordon scalar which forms a 
bound state in the gravitational and electrostatic potentials sourced by its 
own probability current. We begin by specializing the Lagrangians, the field 
equations and the normalization condition to the bound state ansatz f lSOt - 
|52|) . We then note that one of the equations can be solved exactly for the 
negative energy gravitational degree of freedom. Substituting this solution 
into the action gives a functional of the remaining fields whose minimization 
yields the remaining field equations. This functional serves as the basis for a 
variational formulation of the problem. The section closes with a numerical 
solution. 
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5.1 Field Equations 

If a bound state forms we can simplify the fields to take the form (I50ti52p . 

= e-^^*F(r) , A, = -^r)6l, (89) 

goo = -B{r) , fi-oi = , gij = TCij + fif,A(r) . (90) 

Here fj = Xj/r is the radial unit vector and vTjj = 6ij — ViVj is the transverse 
projection operator. The nonvanishing components of the affine connection 
are, 



B' 



B' 



Oi 



A' 



^2B^ ' 

The nonzero components of the Riemann tensor are, 



B" B'^ A'B 



R 



iOj 



jki 



+ 



+ 



2B 4fi2 4AB 
A' 



B' 



''''^2rAB^ 



V rA J 

(91) 



(92) 



TTikrjri-TTurjrk 



^2rA^ ^ 



A' 



+ 



'A-l 



(^) ■ (93) 



Contraction gives the Ricci tensor R^y = R^'^p^, 

/2 



Roo 
R 



B" B" _A'B' A' 
2A~AAB^^A^^Va ' 

B" B'^ A'B' A' 



+ 



+ 



+ 



2B 4^2 AAB rA 



B' A' (A-\ 



2rAB 2rA^ V r^A 
And contracting that gives the Ricci scalar R = g^^R^y, 



R 



B" B" A'B' 2 A' 2B' 



AB 2AB^ 2A^B rA^ rAB 



(94) 



(95) 



(96) 



In specializing the various Lagrangians fll5m7p to our ansatz f lSUll^Ul) it is 
best to use spherical coordinates and integrate the angles so that the effective 
measure factor is, 

= A'n:r^JA{r)B{r) . 



(97) 
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ADM long ago worked out the surface term one must add to the Hilbert 
Lagrangian fllSp for asymptotically flat field configurations [22] . For our 
static, spherically symmetric geometry the result is, 



ADM 



, d (B\ 

-9 ' — 



-9 



(9r V AB 



8nG 



A' 



j.j^2 V j-'^A 



The other two Lagrangians (fTGlfTj 
''-EM > 7rT7TV~9 } 



require no surface terms, 



B 



Hence the action is, 
S[F,<l>,B,A] = {U-t_ 



2AB 



A^ ' 



MlF\ 



-9 ■ 
(98) 

(99) 
(100) 



dr y/—g 
v2 



SttG 



A' 



rA^ 



+ 1 r-^A ' 



2AB 



+ 



B 



F'^ 1 

— -M',F'^.{101) 



It is well known that the operation of making an ansatz for the solution 
does not commute with varying the action to obtain the field equations. 
The correct procedure is to vary first and then simplify. However, all of the 
equations one gets by simplifying first and then varying are correct, and it 
turns out that the only ones we miss are the trivial relations implied by the 
Bianchi identity |23]. Except for the overall factor of At = (t+ — t-), the 
variations with respect to various fields are. 



_1_6S_ 
AtJF 
1 6S 
At 5$ 
_1_6S_ 
AtW 



_1_5S_ 
AtSA 



d_ 
dr 
d_ 

dr 



2F' 



-9 



A 

AB 

" A' 



+ 



2{E + e^fFy/^ 
B 

2e(E + e$)F2y=^ 



2MlF^ 



SnG 



rA^ V f2j\^ 



+ 

{—] 



B 

2AB 
{E+e^fF' 



(102) 
(103) 



n/2 



B 



A 



- M^F^ 



2B 



,(104) 



8nG 



B' 



rAB 



"^V r2 4 ) 



2AB 
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Our goal is to find F{r), $(r), B{r) and A{r) so as to make each of tlie 
variations (110211105]) vanisli. However, tfiose equations fiave solutions for any 
constant E. As always, it is normalizability which puts the "quantum" in 
quantum mechanics. For our problem the normalization condition is, 

dr^g ^ = 1 . (106) 

Any field configurations -Fo('")) '^o('")) Bo{r) and Ao(r) which make all the 
variations f ll02til05]) vanish and also obey (11061) will only do so for very special 
values of E. Note that the field equations and normalizability requires A{r) 
to approach one at infinity, and F{r) to approach zero. The asymptotic 
values of $(r) and B{r) are both gauge choices. By choosing the U{1) gauge 
parameter to be —t x we can make $(r) vanish at infinity. By changing 
time to t/k we induce the rescalings, 

B{r) — > P X B{r) , (107) 
$(r) — > kx <l>(r) , (108) 
E — > kx E . (109) 

We shall always use this freedom to make B{r) approach one at infinity. 

If we can find a normalized solution FQ{r), ^o{r), -Bo(^) and Ao(r) then 
our zeroth order result for the mass is. 

Moth = E- i^^[Fo, $0, Bo, Ao] . (110) 

The first term on the right hand side of (11 101) is from the the two scalar fields, 

0(t+,x)0*(t_,O) Fo(r)Fo(0)e-^^^* . (Ill) 

The final term in (11 101) represents the gravitational and electromagnetic con- 
tribution to the mass. Note that the scalar action vanishes for solutions 
because the scalar Lagrangian is a surface term which goes to zero, 

-d-r —J^o ■ ^ ^ 
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5.2 A Variational Formalism 

Solving differential equations is tough, and we are not able to find exact 
solutions for all four of the fields. For many bound state problems in quantum 
mechanics the absence of exact solutions is not crippling because variational 
techniques allow one to derive strong bounds on the ground state energy. 
Such a technique would be simple to formulate for our Klein-Gordon scalar 
if only the electromagnetic and gravitational potentials were fixed. However, 
the fact that these potentials are sourced by the Klein- Gordon wave function 
itself endows this problem with a slippery, nonlinear character. The presence 
of gravitational interactions is especially problematic because some of the 
constrained degrees of freedom in gravity possess negative energy. Instability 
is only avoided by constraining these degrees of freedom to obey their field 
equations; attempting to minimize the action with respect to these degrees 
of freedom would carry one away from the actual solution. 

There are good reasons for suspecting that the field B{r) is the only 
negative energy degree of freedom. In the normal ADM formalism B = N"^ 
would be the square of the lapse field, and it could be specified arbitrarily 
as a choice of gauge. However, 5 is a dynamical degree of freedom in this 
problem. The structure of our Lagrangian is similar to the usual formalism 
for describing cosmological perturbations during primordial inflation [23] . In 
that setting, as for us, the Lagrangian can be written as the sum of a "kinetic" 
part K and a "potential" part P, 




(113) 



For us the kinetic and potential parts are. 



K 



2A 



(114) 



P 




(115) 



The field equation for B is algebraic and has a trivial solution. 




B 



K 



P 
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Substituting flll6p into f lll3p allows us to express the action in terms of 
just F{r), $(r) and A{r), 



Jo 



(117) 



It is simple to show that varying f lll7p gives equations f ll02lll03p and f llOSp . 
Because (11170 is positive semi-definite, the problem of extremizing it is likely 
to be the same as that of minimizing it. The corresponding normalization 
condition is, 



lAK 



And our zeroth order result for the scalar mass becomes, 
Moth 



E-8ttJ^ drr^^jAoKoPo . 



;il8) 



(119) 



We illustrate the method with simple trial functions for -F(r), $(r) and 
A{r). We cannot make F{r) a spherical shell like ADM, or even a hard sphere, 
because the factors of F'"^ become ill-defined if F{r) has a discontinuity. The 
next best thing is to assume the scalar drops linearly to zero within some 
distance R, 

F{r) = a{R-r)9{R-r) . (120) 
Comparably simple forms for the potentials are. 



$(r) 
1 

I(^ 



ee{R-r) _ ee{r-R) 



e{R-r) + 



c d 
1-- + - 



e{r-R) 



We can make A{r) continuous by choosing, 

c d 
^R^~R^' 

The corresponding forms for the kinetic and potential terms are. 



(121) 
(122) 

(123) 



K{r) = [E-^)\'{R-rmR-r) + £-^{l-^ + ^)e{r~R), (124) 



P(r) 



36 
8^ 



l-br^ + M^{R-rf 



e{R-r) + 



dd{r-R) 



(125) 
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Hence the potential B{r) is 
B(r) 



{E-^fa\R-rfe{R-r) , ^(l-£ + i|)g(r-i?) 



+ 



(126) 



Enforcing that B{r) goes to one at infinity determines the coefficient rf, 

d = aG , (127) 

where a = e^/47reo ~ 1/137 is the fine structure constant. Enforcing conti- 
nuity at r = — which means B{R) = — requires the choice, 



„ aG 



(128) 



At this stage the free parameters in our trial solution are R, a and the 
energy E. The next step is to enforce normalizability, which requires, 



1 

Stt 



{E-lfa\R-rf 



Vl-&r2 J ^_a2^l-hr^+M^{R-ry^ 



{RE-afR'a' 



1 dx 



y/T^ 



X^ 



(129) 
(130) 



-[i-x^+Rm^{i-x)'^ 

= {RE-afR^a^ X l(GR^a^,R^Ml) . (131) 



2\ 



The function /(x, y) in equation (11311) can be expressed in terms of elliptic 
integrals but we may as well treat it as an elementary function and use it to 
express the energy. 



1 



E 

/8nR^aH{GR^a^,Rm^) 

We can now regard the two free parameters as, 

Ra 



+ a 



1 

R 



(132) 



A 



Mn 



/i = RMo . 



The field action is, 



S_ 
At 



/8ttAJ{GM^A^,I2 
JiGMlA^^ 



a 

+ - 



Mo, 



(133) 



(134) 
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where the new integral is, 



J{GM^A\^i 



87rGMoM2 



[l-x2+/i2(l-a;)2] . (135) 



We determine the free parameters A and yU by minimizing S/At. 

Expressions fll3ip . f ll34p and f ll35p seem very complicated. However, 
note that because I{GMqA^,ii^) is an increasing function of A and /i, and 
JiGM^A^ , jj^) is a decreasing function of A and /x, S'/At is a decreasing 
function of /i at fixed A. Hence S/At is minimized, at fixed A, by choosing 
/X to be the maximum value for which the two integrals remain real. 



87rGMoM2 



1 



Q< A<Ar 



At fi = /imax the two integrals become. 



I{GM^A\fji 



''max) 



-X] 



J{GM^A\pi 



'max) 



ov/T^^<,Jl-(l-a:)2]+x2 
i(ixx^(l — x) 



Vl^ 

And the field action (11341) takes the form, 

j (/^max) 



/^max[l-(l-a;) 



(136) 

(137) 
-x2 (138) 

(139) 

^l + /^max)«(^n 

At this stage the problem is numerical. Fig.[5]shows S/M^At as a function 
of /Umax for GMq = 0.36. The minimum seems to be at about /imax = 0.12, 
which corresponds to, 

A'max 



At 



a 



GMi 



Mn . 



R 



GMi 
l3GM^ 



xVG^ 0.20 X VG 



pmaxA/ J- 
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We will see that these results are not very accurate. 
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Figure 5: Plot of j-^^ versus //max from equation fll39p for GM^ = 0.36. 
The minimum seems to be at about /imax = 0.12. 

5.3 Numerical Results 

It is desirable to check any variational ansatz against a direct, numerical 
solution to the problem. Of course computers can only solve for dimensionless 
quantities, so it is first necessary to express everything in geometrodynamical 
units, using G to absorb each quantity's natural units, 

— MEM 
Gf,Mo = -=,i? = -=,M = -=, (144) 



F{r) = ^ , $(r) = M!:)^ , Bir) = B{r) , A{t) = A{r) 
\/G e\/G 



(145) 

v/^ = g/^, K{r) = ^ , P{r) = ^ . (146) 

Note that we have absorbed the energy into the electrostatic potential. In all 
cases we employ a tilde to denote the dimensionless quantity. Geometrody- 
namic fields such as F are considered to be functions of the geometrodynamic 
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Figure 6: Plots of the scalar amplitude F{r) (in units of Mpi = l/y/G) and 
the electrostatic potential $(r) (in units of Mpi/e) as functions of r (in units 
of 1/Mpi). These figures were generated for bare mass Mq = 0.60 Mpi. 



radius r. A prime on such a field indicates differentiation with respect to r, 
so we have, 



F' $' 
F' = — $' = — 

G ' eC 

In these units the four field equations fll02til05|) take the form, 
+ - MlFJ-g = , 
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(Recall that a = e'^/AireQ ^ 1/137 is the fine structure constant.) The kinetic 
and potential terms are. 
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Figure 7: Plots of minus the tt component of the metric B{r) (dimensionless) 
as a function of r (in units of 1/Mpi = a/G). The right hand figure has an 
expanded vertical axis to show the small variation of the field. These figures 
were generated for bare mass Mq = 0.60 Mpi. 
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The normalization condition is, 
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And the final result is, 
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The nonlinear nature of this problem requires a special solution strategy. 
The development of our technique was facilitated by the vast amount of work 
that has been done of "boson stars" [251 126] . There has also been a recent 
study by Carlip of gravitationally bound atoms [27] . 

Our strategy is to begin by evolving equations fl 148111 51]) outward from 
f = 0, with arbitrary choices for F{0) > 0, $(0) < and -B(O) > 0, and with 
the other boundary values at, 



F'(0) = , $(0) = , A{0) = 1 



(156) 
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Figure 8: Plots of the rr component of the metric A{r) (dimensionless) as 
a function of r (in units of 1/Mpi = y/G). The right hand figure has an 
expanded vertical axis to show the small variation of the field. These figures 
were generated for bare mass Mq = 0.60 Mpi. 

The choice of -B(O) > really is arbitrary because we will eventually make a 
global re-scaling of time to force B{r) to approach one at infinity. However, 
the choice of $(0) essentially gives the energy, and this matters of course. 
There is zero probability of guessing a true eigenvalue. With the other con- 
ditions fixed, varying $(0) gives solutions for which F{r) either becomes 
negative (which a magnitude cannot do) or grows at infinity (which a nor- 
malizable solution cannot do). One knows that a true energy eigenvalue has 
been bracketed between two different choices of $(0) when the behavior of 
F{r) changes from one extreme to the other. Then one closes in on the eigen- 
value to whatever accuracy is desired. Note that this means cutting off the 
behavior of the solution past a certain value of r, beyond which F[r) begins 
to degenerate. 

The procedure we have just outlined gives a solution which is normal- 
izable, but not yet normalized. For that we compute fll54p and then either 
increase or decrease F{0) as needed. Of course the nonlinear nature of this 
problem means that one does not get a solution by simply multiplying F{r) 
by a constant! We must instead start from the new -F(O) and again go 
through the process of trapping the energy eigenvalue. However, our evolu- 
tion programs are efficient enough that this can be done to high accuracy. 
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and fairly quickly. 

Figures EHHl show the behavior of the fields for M = 0.6. For this bare 
mass the energy is E ^ 0.586378 and the total mass is M ^ —0.000882. A 
measure of the numerical error is the accuracy with which the scalar action 
vanishes, which is 5'sc/At 3 x 10~^. Another measure of accuracy comes 
from the finite cutoff at r = Rent, occasioned by the finite accuracy of E. 
For that bare mass we cut the various integrations off at i?cut = 500, which 
corresponds to a contribution of a /Rent ~ 10~^ from the electromagnetic 
tail. 

The variational results fll40H143p obtained in the previous subsection are 
quite different from the numerical solution. From Fig. one can see than 
the scalar amplitude has roughly the same shape as that of our trial function 
(11201) . but with initial height Ga ~ 0.0065, rather than the variational value 
f|T4Tl) of Ga ^ 1.7. And the radial extent is about VG R ^ 300, rather 
than the variational result f ll40p of a/G R ~ 0.20. There is no chance that 
the discrepancy derives from the numerical solution, the error of which we 
estimate to be no larger than 10~^. The problem must lie instead with the 
variational formalism. Our trial solution seems roughly correct, but it may be 
that, like B{r), the gravitational potential A{r) represents a negative energy 
direction in field space. In that case minimizing the constrained action would 
take us away from the actual solution, which seems to be what has happened. 

Table [Ugives our results for E and M for a variety of different bare masses. 
The most obvious feature is the almost total cancellation between the energy 
of the scalar wave function and the field action, to give a very small, negative 
total mass. This is physical nonsense because it fails to agree with the mass 
one can read off from asymptotic values of the metric. We believe that the 
problem arises from the asymptotic conditions f l?ni^2|) — which are certainly 
valid for scattering with other particles — not being right for the study of 
self-interactions. We believe that this can be fixed without much change. 

The other features of our numerical work are: 

• The energy E agrees with the mass inferred from the asymptotic values 
of the metric. 

• There is no bound state unless the bare mass Mq exceeds the ADM 
result of ^/a ^ 0.85 [26] . 

• The bound state energy E is in all cases less than the bare mass. 
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M 






0.15 


0.149994 


-0.000320 


2.0 X 10"^° 


70,000 


0.20 


0.199962 


-0.000260 


1.5 X 10"^ 


25,000 


0.25 


0.249870 


-0.000505 


2.9 X lO-'' 


11,750 


0.30 


0.299653 


-0.000239 


9.6 X lO-'' 


4,500 


0.35 


0.319215 


-0.001075 


5.7 X 10-'^ 


3.200 


0.40 


0.398424 


-0.000455 


1.3 X 10"^ 


1,800 


0.45 


0.447073 


-0.000826 


3.5 X 10-^ 


1,300 


0.50 


0.494904 


-0.000788 


7.3 X 10-^ 


1,400 


0.55 


0.541546 


-0.000412 


1.6 X 10-6 


1,100 


0.60 


0.586378 


-0.000815 


3.4 X 10-6 


500 


0.65 


0.628511 


-0.001440 


8.1 X 10-6 


350 


0.70 


0.666426 


-0.001866 


1.7 X 10-5 


275 


0.75 


0.696992 


-0.000532 


4.1 X 10-5 


400 



Table 1: Numerical results for tlie^ scalar energy E and the total mass M for 
different values of the bare mass Mq. Also given are the scalar action S'sc/Ai, 
which should vanish, and the cutoff radius beyond which the finite accuracy 
of the energy eigenvalue makes the solution unreliable. All quantities are 
expressed in Planck units. 
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• The ratio E/Mq increases with Mq and eventually becomes zero |26] . 



Because there would not even be any bound states without gravity, it seems 
fair to conclude that the system depends nonanalytically upon G. 



We have explored the possibility that the apparent problems of quantum 
general relativity may be artifacts of conventional perturbation theory. One 
might think this unlikely because the absence of recognizable, low energy 
quantum gravitational phenomena implies that some asymptotic series ex- 
pansion is wonderfully accurate. However, it may be that the correct series 
involves logarithms or fractional powers of Newton's constant. If that were 
the case, trying to re-expand in integer powers of G would result in an esca- 
lating series of divergences, which is exactly what conventional perturbation 
theory shows. 

We studied this possibihty in the context of computing the mass of a 
charged, gravitating scalar. An exact result for the classical limit of this 
system was derived by ADM in 1960 [15], and it does exhibit both nonanalytic 
dependence upon G and the breakdown of conventional perturbation theory. 
If the classical point particle is regulated to be a spherical shell of radius R, 
the ADM result is. 



The correct zero radius limit is M = ya/G. Its finiteness results from neg- 
ative gravitational interaction energy canceling the positive electromagnetic 
energy. In contrast, the perturbative result is obtained by first expanding 
the square root in powers of G and e^, which produces a series of ever-higher 
divergences with alternating signs. The alternating signs are a signal that 
gravity is trying to cancel the electromagnetic self-energy divergence, but this 
cancellation can never happen in conventional perturbation theory because 
the gravitational response to a divergence at one order is delayed until one 
order higher. What we need for quantum gravity is an alternate expansion 
in which the negative gravitational interaction energy has a chance to "keep 
up" with what is going on in the positive energy sectors. 



6 Epilogue 




(157) 
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In section 3 we derived an exact functional integral expression fl4Up for the 
scalar mass. We then developed an alternate asymptotic expansion based on 
the Method of Stationary Phase, with the full functional integrand — not 
just the action — used to determine the stationary point. This is more diffi- 
cult to implement than conventional perturbation theory, but it is also more 
correct. A simple integral representation for the Bessel function illustrates 
the distinction between our approach and that of conventional perturbation 
theory, 

J^[z) = ^ rrf0e*^^*°(^) X (e-*')"^ . (158) 



27r 

In our approach both factors are included in the exponent and the two sta- 
tionary points are found by minimizing the function f{9) = zsin(6') — A^^^, 

f\e±) =0 =^ 9± = ±acos(— ) . (159) 
The values of the function and its second derivative at these points are, 

/(^±) = ±[Vz^-N^ - Aracos(^)] , f"ie±) = ^Vz^-N^ . (160) 
And the result for the 0th and 1st order contributions is. 



Jn(z) — > J — , ^ cos\Vz^-N^ - Nacos(—)--] . (161) 



In contrast, conventional perturbation theory would be based on the function 
f{9) = zsm{9), with the stationary points at 9± = ±|. The result for the 
0th and 1st order contributions from conventional perturbation theory is. 
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TTZ 



TT TT 



Jn{z) > W — COS Z-N . (162) 



2 4 



Section 4 presents an analysis of the new expansion in the context of 
a simplified model. We conclude that all the old £ loop diagrams appear 
at £-th order in the new expansion. However, the old i loop diagrams are 
combined with an infinite class of new diagrams which possess more external 
lines and no more than i loops. The new diagrams which are added at i-th 
order are all subtracted at higher orders, so we are really adding zero to the 
usual expansion. Because the new ^-th order diagrams have no more than 
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(. loops, the divergences of the new expansion can be no worse than those 
of conventional perturbation theory. Because infinitely many new diagrams 
are added at each order, the new expansion can depend nonanalytically on 
Newton's constant. It also offers a way in which the negative gravitational 
interaction energy can respond, at the same order, to problems in the positive 
energy sectors. These are all desirable features, although it must be admitted 
that these is no guarantee at this stage that the new expansion is any better 
than the old one. 

The analysis of section 4 was done only to understand how the new ex- 
pansion compares with the old one. There are much better ways of actually 
implementing the new expansion. We exploit two of these methods in section 
5 to evaluate the zeroth order result. Our analysis is based on interpreting the 
zeroth order term as the phase developed by a first-quantized Klein-Gordon 
scalar moving in the gravitational and electrodynamic potentials which are 
sourced by its own probability current. The fact that this system reduces 
to the ADM problem for — )■ provides a solid reason for believing both 
that the negative energy gravitational interactions cancel at least some of 
the usual self-energy divergences, and that the final result depends nonana- 
lytically on Newton's constant. 

Evaluating the zeroth order term of the new expansion amounts to solving 
for a bound state of the scalar in its own potentials. Although we cannot 
obtain exact solutions for all four of the relevant field equations, we were 
able to eliminate one of the negative energy gravitational degrees of freedom 
to derive a variational formalism. We were also able to solve the equations 
numerically, taking advantage of the vast body of work which has been done 
on "boson stars" [25|[26]. 

We achieved high numerical accuracy which revealed a substantial dis- 
crepancy with the variational approach. This probably means that the grav- 
itational potential we were not able to eliminate also carries negative energy, 
so that minimizing the constrained action takes one away from the true so- 
lution. As was seen in previous numerical work, we found that there are no 
solutions unless the bare mass Mq is greater than the ADM result of \J a/G. 
We developed solutions for many choices of Mq above this limit. All of them 
show an almost total cancellation between the energy of the scalar wave func- 
tion and the field energy of the gravitational and electromagnetic potentials. 
This gives nearly zero for the total mass, which seems to be nonsense. It also 
fails to agree with a determination of the scalar mass from the asymptotic 
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values of the gravitational potentials. 

The problem seems to derive from our use of the asymptotic conditions 
f l2ni22|) . Expressed in simple words, these conditions mean that "the fields 
become free at asymptotically early and late times". That is perfectly true 
(in the weak operator sense and assuming the existence of a mass gap) for 
interactions between different particles, which is the usual application |TS] . 
However, we are here trying to use the conditions to study interactions of a 
particle with itself. These self-interactions would usually be subsumed into 
forcing the field strength and mass to come out right by renormalization, but 
that is exactly what we are not doing. We believe that when a more accurate 
procedure is used to interpolate the single particle states — which might be 
as simple as including a U{1) gauge string between the two fields to make 
them invariant — then the nonsense result for M will go away, and most of 
our analysis of the new expansion will be unchanged. 
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